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Abstract 

Analytical  results  for  complex  moduli  of  uniaxially  fiber  reinforced  mat¬ 
erials  made  of  viscoelastic  matrix  and  elastic  fibers  are  reviewed.  A  general 
method  is  established  to  predict  complex  moduli  and  loss  tangents  of  viscoelas¬ 
tic  laminates  made  of  uniaxially  reinforced  laminae.  Application  of  results 
is  demonstrated  by  two  examples  of  analysis  of  damping  of  structural  vibrations: 
Attenuation  of  vibrations  in  uniaxially  reinforced  Timoshenko  beam  and  torsional 
vibrations  of  laminated  cylinder. 


^  Presented  at  8th  International  Council  of  the  Aeronautical  Sciences  (ICAS) 
Congress,  Amsterdam,  August  28  -  September  2,  1972. 

Professor..  Department  of  Materials  Engineering, 

Technion  -  Israel  Institute  of  Technology. 
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1.  Introduction 


The  ever  increasing  use  of  fiber  composites  for  aero/space  structures 
requires  the  development  of  rational  methods  for  prediction  of  their  rele¬ 
vant  properties  within  engineering  accuracy. 

In  most  current  fiber  composites  the  matrix  is  a  polymer  such  as  epoxy. 
It  is  well  known  that  such  polymers  exhibit  the  effect  of  vibration  damping. 
Therefore  such  damping  effects  will  also  occur  in  composites  in  which  the 
matrix  is  polymeric.  Since  aero/space  structures  are  subjected  to  severe 
vibrational  environment  and  since  vibration  damping  is  beneficial,  the 
quantitative  prediction  of  such  damping  is  of  considerable  engineering  imp¬ 
ortance. 

It  is  of  interest  to  emphasize  the  unique  combination  of  desirable  pro¬ 
perties  which  are  exhibited  by  fiber  composites:  Superior  strength  and 
stiffness,  low  weight  and  vibration  damping.  No  other  materials  seem  to 
possess  this  many  advantages. 

In  order  to  handle  the  problem  analytically,  it  is  assumed  that  the 
matrix  is  linearly  viscoelastic.  Its  dynamic  viscoelastic  properties  can 
then  be  characterized  in  terms  of  the  usual  complex  moduli  of  viscoelasticity 
which  are  assumed  to  be  known  on  the  basis  of  experiments.  The  fibers  are 
represented  as  linear  elastic.  The  composite  with  such  constituents  behaves 
macroscopically  as  a  linear  viscoelastic  body  which  is  characterized  by 
effective  complex  moduli.  There  arise  three  classes  of  important  investi¬ 
gation: 
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(a)  Prediction  of  effective  complex  moduli  of  a  uniaxial ly  reinforced 
material  on  the  basis  of  matrix  complex  moduli;  fiber  elastic  moduli  and  in¬ 
ternal  geometrical  parameters  such  as  constituent  volume  fractions,  fiber 
shapes,  etc. 

(b)  Prediction  of  the  effective  complex  moduli  of  a  laminate,  whose 
laminae  are  composed  of  uniaxial ly  reinforced  material,  on  the  basis  of  the 
uniaxial  material  effective  complex  moduli  found  in  (a)  and  the  laminate 
internal  geometry. 

(c)  Viscoelastic  vibration  analysis  of  structures  made  of  fiber  compo¬ 
sites. 

These  different  kinds  of  problems  will  be  discussed  consecutively. 
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2.  Complex  Moduli  of  Uniaxial ly  Fiber  Reinforced  Materials 

A  general  theory  of  prediction  of  effective  complex  moduli  of  composites 
with  linear  viscoelf  Stic  constituents  has  been  given  previously  [1],  [2],  [3], 
It  will  here  suffice  to  discuss  without  proof  some  results  which  are  pertinent 
for  the  present  investigation. 

Let  the  local  average  strains  and  stresses  in  a  composite  be  of  oscillatory 
nature.  Thus: 


L*  >  *”■  £••  6 

13  13 


not 


(2.1) 


a. .  =  a. .  e 
13  i3 


Udt 


where  overbars  denote  average,  l  =  /T,  «  is  frequency,  t  is  time  and  latin 
subscripts  range  over  1,  2,  3.  The  effective  complex  moduli  C^jj,  of  a  generally 
anisotropic  composite  are  defined  by  the  relation: 


°ij  a  ^ijk£(uo)ek2. 


'  '  clSki« 


(a) 


(2.2) 


(b) 


where  superscripts  R  and  I  denote  real  and  imaginary  parts  respectively. 

The  assumption  is  made  that  the  fiber  reinforced  material  under  considera¬ 
tion  is  macroscopically  transversely  isotropic  with  respect  to  fiber  direction. 
Then  (2.2a)  assumes  the  form: 


*• -^s^i5BSiSSH^«y^^»&6a^»»«PSeras*!^«S*SM^?@aS!3SBS®S»»^g5SgB6igSSH^S#S59P3fStR*«3a 
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0u  =  Cu  £u  ■*■  C 1 2  C22  +  C12  £33 


022  =  C12  Ell  +  C22  £22  +  C22  £32 


033  s  C12  Eli  ♦  C23  £22  +  C22  £3- 


012  =  2C44  &12 


(2.3) 


023  *  (C22  "  C23)£23 


03 1  =  2C44  £31 


where  xi  is  in  fiber  direction  and  X2,  X3  are  in  the  transverse  plane.  Fig.  1. 


In  another  notation,  the  complex  moduli  in  (2.3)  are  written: 


Cn  *  n 


C12  =  A 


C22  -  k  ♦  G» 


(2.4) 


C23  =  k  -  G- 


C44  =  G, 


Here  k  is  a  transverse  complex  bulk  modulus,  G,p  -  transverse  complex  shear 
modulus  in  X2,  X3  plane  and  G^  -  axial  complex  shear  modulus  in  planes  contain¬ 


ing  fiber  direction  xi.  The  physical  interpretation  of  n  and  l  is  here  of 


little  interest. 


Inversion  of  (2.3)  is  written  in  the  form: 


■*-  «frWf*  -W/SV t,  r£  <  ■ 


P  6 


6  - 


1  VA  VA 

Ell  *  ~~  3u  -  ~  022  ~  ~  033 

ea  ea  E  a. 


~  A  ~  1  ~  T  ~ 

e223-  —  on  +  — -  022  -  fr*°33 

“A  eT  eT 


v  \j 

~  A  -  T  ~  1  ~ 

C33=-  ~~0n  -  xr-  022  *  ~  O33 

ea  et  Ex 


£12 


E2  3 


£31  a 


012 

2GA 


O23 

2Gx 


03: 

22a 


(2.5) 


where  E^  and  F^.  are  complex  Youngs’  moduli  in  axial  (fiber)  direction  and 
transverse  (to  fiber)  direction,  respectively,  and  vA  and  Uj.  are  associated 
complex  Poisson's  ratios. 

Establishment  of  analytical  expressions  for  the  various  effective 
complex  moduli  listed  above  in  terms  of  matrix  complex  moduli,  fiber  elastic 
moduli  and  phase  geometry  is  based  on  a  correspondence  principle  [1] ,  [2] 
which  states:  The  effective  complex  modi  1  of  a  viscoelastic  composite  are 
obtained  by  replacement  of  phase  elastic  moduli  by  phase  complex  moduli 
of  a  composite  with  identical  phase  geometry. 

In  the  usual  fiber  reinforced  materials  the  following  conditions  are 
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usually  fulfilled  with  sufficient  accuracy 

(a)  Fibers  are  by  an  order  of  magnitude  stiffer  than  matrix 

(b)  The  matrix  is  isotropic  and  is  viscoelastic  in  shear  only 

(c)  The  matrix  shear  loss  tangent  is  not  larger  than  0.1.  Thus: 


tan  6"1  =  *  0.1 

Gm 


(2.6) 


Under  these  conditions  the  following  results  have  been  shown  [3] 
to  be  valid. 

(d)  The  imaginary  parts  of  the  effective  complex  moduli,  n,  1,  ic, 
Ea  are  much  smaller  than  the  imaginary  parts  of  the  effective  shear 
moduli  Ga,  Gt,  and  of  E^.. 

(e)  To  obtain  real  parts  of  all  effective  complex  moduli  it  is 
merely  necessary  to  take  corresponding  expressions  for  effective  elastic 
elastic  moduli  and  to  replace  in  them  matrix  elastic  moduli  by  real 
parts  of  matrix  complex  moduli. 

Some  simple  general  results  which  are  valid  under  the  conditions 
listed  above  will  now  be  give*: 

eaC«0  -  eb(h»)  .  Ef  vf 


E>)  -E>)vB.Efvf 


bJCu)  ■  B‘(.)  V, 

.1 


m 


(2.7) 


tan  6_  =  — =- 
E  _R 


tand 


m 


1  ♦ 


WE>, 


«  tan  6, 


m 


m 


y,/*iwr*ir 
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Here  a  and  f  denote  mati'ix  and  fibers,  respectively,  and  v  stands  for 
volume  fraction,  tan  6£  is  the  loss  tangent  for  uniaxial  stressing  in  fiber 
direction  while  tan  6®  is  the  corresponding  loss  tangent  for  the  isotropic 
matrix. 


For  axial  shear 


1  +  v# 


GA(tw)  =  Gm(ito)  y 


ej(w) 


gJ(«) 


Jt  1  4  Vf 

Ja  1  -  v. 


4  1  +  vf 
“n  1  ■  v. 


(2.8) 


tan  6, 


tan  5 


m 


For  transverse  shear 


G*  =  g!J  tan  6® 


tan  6„  *  tan  5' 


m 


(2.9) 


,R 


An  expression  for  G^,  has  been  given  in  [2] ,  [3] .  Results  for  other 
complex  moduli  may  also  be  found  in  these  references. 


5?SsrtRHS^ 

! 


SKS^^^’S^'R’iSv'' 


3.  Complex  Moduli  of  Laminates 


The  laminates  to  be  considered  are  composed  of  plane  laminae  of 

*:  ) 

uniaxially  fiber  reinforced  materials ■,  the  direction  of  reinforcement 

:  1  5 

being  different  in  each  lamina.  The  laminate  is  referred  to  a  fixed 
coordinate  system  x»,  X2,  xs  where  xi,X2  are  in  the  plane  of  the  lami- 

i 

nae  and  X3  is  normal  to  it,  Fig.  2.  The  nth  lamina  in  the  laminate  is 
referred  to  a  material  system  of  axes  ,xi^,  X2^,  xj  where  xi^  is 
normal  to  the  fibers  and  X3  coincides  with  the  laminate  X3.  The  posi- 

l  , 

tion  of  the  xx^n\  X2in^  system  is  defined  with  respect’ to  the  xi,  X2 


system  by  the  reinforcement  angle 


en  *  *  (xx(n),  Xj) 


(3.1) 


Fundamental  assumptions  of  fiber  composite  laminate  theory1  are:  (a)  Any 

1 

lamina  can  be  replaced  by  a  homogenous  material  whosei  properties  are  the 
effective  properties  of  the  uniaxial  FRM  o£  which  the  lamina  is  made. 

(b)  The  laminae  are  in’  statesof  plane  stress. 

First  an  elastic  laminate  will  be' considered.  The  plane  stress-strain 
relations  of  a  lamina  referred  to  its  material  system  of  axes  xi^,  X2^n' 


are  then: 


„  (n)  1  (n)  VA  _  (n) 

Ell  “  g-Ou  -  022 

CA  A 


0  (n)  VA  (n)  1  „  (n) 

G22  *-  r-ffll  ♦  p—  O22 


(n)  „  012 


(3.2) 


**WW«ty 


^^rsissmawuw^^ 
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where: 

E^  -  axial  Youngs  modulus  (in  fiber  direction) 


\>A  -  associated  axial  Poisson's  ratio 


Et  -  transverse  Youngs  modulus  (normal  to  fibers) 


-  axial  shear  modulus  (in  Xj^,  X2^  plane) 


The  inverse  of  (3.2)  is: 


_  (n)  -  _  (n)  -  .  (n) 

Oil  *  Ctl  env  ♦  Cl 2  622 


022  ^  -  C12  6llCn)  +  C22  622 (Tl)  (3.3) 


Oi2(n^  »  2C|»I*  6j 2 


where 


C11  *  Cuu 


C12  *  Ci 122  = 


vaet 


1  -  A  V 

E"  A 
“T 


C22  *  C2222  * 


E»  ..2 


1  -  A  V 


*F 


Chh  a  C1212  *  G. 


(3.4) 


! 

1 

.1 


8 

I 

I 


4 

1 

V, 


'is 

% 


In  terns  of  the  four  index  moduli  in  (3.4),  (3.3)  can  be  written 


compactly  as: 


,(n) 

J  <xBy$ 


(3.5) 


where  here  and  from  now  on  Greek  indices  range  over  1,  2. 

For  the  sake  of  simplicity  there  will  be  considered  the  special 
group  of  laminates  in  which  the  application  of  membrane  force  NQg  in  the 
plane  of  the  laminate  does  not  induce  bending  or  torsion  in  the  laminae. 
The  most  important  kind  of  laminate  which  fulfills  this  requirement  is  a 
symmetric  laminate.  Such  a  laminate  has  the  property  that  its  middle 
plane  is  a  plane  of  symmetry  for  the  geometry  and  elastic  moduli  of  the 
laminate.  The  laminate  is  thus  composed  of  laminae  pairs  in  each  of 
which  the  laminae  are  of  same  thickness,  are  symmetrically  located  with 
respect  to  the  middle  plane  and  have  the  same  elastic  properties  with 
respect  to  the  Xi,  X2  system. 

The  last  condition  is  most  commonly  fulfilled  by  laminae  made  of 
identical  material  and  same  reinforcement  angle  0n  (3.1),  in  each  pair. 
The  elastic  stress-strain  relation  of  such  a  laminate  is  given  by 


°a  rx’  Ca;*SeYS  W 


N 


Ca$Y«  ^^a&YS  Vh  (b) 


(3.6) 


ivy**  2. * 


where 


C* 

La8y6 


-  applied  average  plane  stress 


average  strain 


-  laminate  thickness 


-  effective  elastic  moduli  of  laminate 


-  thickness  of  n  lamina 


-  number  of  laminae 


-  the  laminae  elastic  moduli  (3.4)  transformed  to  the  xi,  xz 
system  of  axes. 


For  establishment  of  the  result  (3.6)  see  e.g.  (4J.  Proof  that  (3.6) 
is  based  on  an  elasticity  solution  which  is  exact  in  the  Saint  Venant  sense 
for  a  laminate  whose  thickness  is  small  compared  to  its  plane  dimensions 
has  been  given  oy  B.W.  Rosen  (unpublished). 

Let  it  be  now  assumed  that  the  laminate  is  viscoelastic  but  remains 
symmetric  as  described  above.  The  laminate  is  subjected  to  oscillatory 
membrane  loads 


n«b  a  n«3  6 


(3.7) 


The  average  stresses  associated  with  (3.7)  are  then: 

N  N 

-  aB  aB  lwt  %  \ 

aaB  “  T*  3  TT  e  s  aaB  e 


(3.8) 


The  strain  response  of  the  laminate  is 
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eaB  '  eaB  e 


last 


(3.9) 


The  relation  between  5Qg  and  eQg  is  written: 


°a6  -  Cc.BY«  V 


(3.10) 


where  are  the  effective  complex  moduli  of  the  laminate. 


It  follows  by  the  general  correspondence  principle  of  [1]  which  was 


quoted  above  that  C*^  can  be  expressed  in  form  (3.6b).  Thus: 


N 


'w  -  *  (n)c~*,*  (>“)  »■/>> 


n*l 


aByS 


n' 


(3.11) 


The  single  laminae  complex  moduli  tn)Ca6YS  in  (3.11)  are  interpreted  as 
follows:  In  the  material  axes  xj ^ ,  X2^  of  the  nth  lamina  the  stresses 
and  strains  are: 


a  (n)  „  5  (n)  lwt 

°aB  aaB  ® 


(3.12) 


_  (n)  ~  (n)  xwt 

eaB  eaB 


The  relation  between  and  eog^  in  (3.12)  is  of  type  (3.3-3t4) 


Thus: 


l£9rti»9»J'sW  ..r.  .*•  ' 


I 


1 


1 

i 


fl  W.gCn)  g  Cn) 

aS  aSyS  Y« 


(3.13) 


Clin  ®  Cl J 


&r  , 

‘-E* 


C1122  *  C12 


Vr 

et  , 

1  -  J-vf 
ea  A 


(3.14) 


-2222  ■  U22 


1  -!lv? 

ea  a 


Cl 2 12  *  Cm,  ■  G^ 

where  ea*  ^t*  *a  and  G^  are  the  effective  complex  moduli  of  the  uniaxial 
material  which  were  discussed  in  par.  2. 

If  the  complex  stress  strain  relation  (3.13)  is  transformed  to  the 
laminate  axes  xi,  X2  it  assumes  the  form: 


(n)3  .  -  ta)c  .  .  (n)€  . 
aB  aBY$  y$ 


(3.15) 


This  defines  the  v njCR  .  in  (3.11). 
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By  tensor  transformation: 

^Cim  *  (n)Cn  ■  Cucos%0  +  C22sinl,0  + 

n  n 

+  2Ci2Cos20  sin20  +  4C.„,cos20  sin20 
11  n  n  n 

^Ch22  *  ^612  *  (Cji  +  622) cos 2 0  sin2 0  ♦ 

n  n 

♦  Ci2(cos40  +  sin*1©  )  -  4Ci,%cos20  sin20 
H  n  n  n 

^C2222  *  ^C22  *  Cxi  sin  **0  +  C22cosl,0  + 

n  n 

+  2Cx2cos20nsin20n  +  4Cin,cos20nsin20n 

(3.16) 

^Ciii2  »  ^n^Cxi»  »  _  Cxicos30  sin  0  + 

n  n 

♦  C22cos0  sin30  +  Cx 2 (cos 30  sin0  -  cos0  sin30  )  + 

**  ii  n  n  n' 

1 

+  261,1,  (cos 30 sin0  -  cose  Sin30  ) 
n  n  n  n 

(n)  "*  (nl  •. 

C22 1 2*  ‘’C2H  *  -  Cx  xcos0_sin30  + 

n  n 

+  C22COS 30  sin0  +  C x 2 (cos0  sin 30  -  cos30  sin0  )  + 

**  n  n  n  n  n 

♦  2C^i,(cos0nsin30n  -  cos30nsin0n) 
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^Ci212  a  ^Ci,i»  a  (Clt  +  C22)cos29nsin20n 
-  2Ci2COS20nsin20n  +  C*i*(cos20n  -  sin20n)2 

The  preceding  developments  together  with  the  results  for  complex  modu¬ 
li  of  uniaxially  fiber  reinforced  materials  define  the  computation  methods 
of  the  effective  complex  moduli  of  symmetric  laminates  as  expressed  by  (3.11).' 

For  practical  purposes  it  is  frequently  necessary  to  compute  the  effec¬ 
tive  complex  compliances  which  are  defined  as  the  inverse  of  (3.11).  In 
matrix  notation: 

C*  •  S*  »  J  (3.17) 

•v 

where  S*  denotes  the  effective  complex  compliance  matrix  and  J  is  the  unit 
matrix.  Separation  of  (3.17)  into  real  and  imaginary  parts  yields 

C*R  •  S*R  .-  C*1  •  S*1  =  J  (a) 

(3.18) 

C*R  •  S*1*  C*1  •  S*R  =  0 

Great  facilitation  is  achieved  if  it  is  noted  that  in  view  of  (2.6)  the 
second  term  in  the  left  side  of  (3.18a).  can  be  neglected  with  respect  to 
the  first.  It  then  follows: 

C*R  •  S*R  »  J  (a) 

(3.19) 

T  O  R 

s*  =- S*  *  C  •  S*  (b) 


(3.16 

cont'd) 
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Thus,  once  C*  and  C*  have  been  computed  from  (3.11),  (3.19)  define  the 


effective  complex  compliance  matrix  by  simple  real  matrix  operations. 
Another  important  simplification  is  obtained  if  in  the  laminate 


to  each  pair  with  reinforcement  angle  9n  and  thickness  tn  corresponds 


another  pair  with  reinforcement  -0„  and  same  thickness  t  .  It  is  then 
r  n  n 

easily  realized  by  the  form  of  (3.16)  that  all  contributions  to  (3.11) 


of  terms  with  odd  powers  of  cos0n  and  sin0^  cancel  mutually.  Thus  in 


this  event  the  effective  complex  moduli  matrix  (3.11)  has  the  form 


Cllll  Cl 122  0  . 

Cii  C$2  0 

c*  ■ 

C$ 122  C$222  0 

8 

C$2  C$2  0 

0  0  C$212 

o  o  cU 

(3.20) 


and  so  the  laminate  is  macroscopically  orthotropic.  The  situation  just 
described  is  of  frequent  practical  occurrence.  For  example  a  symmetric 


laminate  with  laminae  reinforcement  in  8^  «=  0,90°,  ±  45°  directions. 


1 

1 

$ 

V 

t 
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4.  Structural  Applications 

4.1  Free  flexural  vibrations  of  a 
~~  fiber  reinforced  beam 


As  a  first  example  there  is  considered  the  case  of  free  flexural  vibrations 
of  a  simply  supported  beam  which  is  uniaxial ly  reinforced  in  beam  axis  direction. 
The  purpose  of  the  investigation  is  to  compare  vibration  damping  due  to  matrix 
viscoelasticity  on  the  basis  of  the  usual  theory  which  neglects  the  effect  of 
shear  and  on  the  basis  of  the  more  refined  Timoshenko  theory  which  takes  into  ac¬ 
count  shear  as  well  as  rotatory  inertia.  For  isotropic  materials,  in  which  the 
complex  Young's  modulus  loss  tangent  and  the  complex  shear  modulus  loss  tangent 
are  of  same  order,  the  added  effect  of  shear  and  rotatory  inertia  is  small  for 
vibration  modes  of  low  order  and  for  long  beams.  In  the  present  case,  however, 
where  the  axial  Young's  modulus  loss  tangent  is  by  an  order  of  magnitude  smaller 
than  that  of  the  axial  shear  modulus  (section  2)  the  situation  is  quite  different 
as  will  be  shown  below: 

Considering  only  the  effect  of  flexure  the  differential  equation  of  the 
freely  vibrating  beam  is: 


$ 

X 


2  92w  .  32w 


w  n  w  n  a 

3^+  W*”  0 


V 


w 


(a) 

00 


(4.1) 


where: 


-  complex  axial  Young's  modulus 

I  -  moment  cf  inertia 
A  -  area  of  cross  section 
p  -  density 

Boundary  conditions  of  free  support  are: 

i 

W,  a  O  X  *  o,i,  (4.2) 

Conventional  viscoelastic  vibrations  analysis  shows  that  the  inodes 
of  vibration  are  given  by 


Equ.  (4.4b)  is  a  repetition  of  (2.7d).  The  results  are  valid  for  small 

enough  loss  tangents,  of  order  (2.6).  The  attenuation  p  is  defined  by 

R 
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Next  the  same  beam  is  considered  in  Timoshenko  fashion,  with  shear 
and  rotatory  inertia.  By  the  correspondence  principle  for  viscoelastic 
vibrations  [5],  the'  elastic  Timoshenko  beam  equation  [6]  transforms  into 


'  cix^  W? 


kE. 

gA 


94w 


kr2 


*  0  (4.6) 


where  c2  is  given  by  (4.1b),  G^  is  the  axial  complex  shear  modulus  (2.8a), 
k  is  the  strength  of  materials  shear  shape  factor  of  the  section  and 
r2  ■  I/A. 

Equ.  (4.6)  with  boundary  conditions  of  simple  support  admits  a  solution 
of  the  form,  [6]: 


w«(x»w  “  A  sin  — -  e  n 
n  n  i 


where  ©n  is  a  solution  of  the  complex  frequency  equation: 


kE. 


rnr 


n  % 


The  solution  of  (4.8)  is  the  complex  "frequency" 


fii  »  Cp'  + 
n  n  n 


(4.7) 


c2a  4-  [1+a  2r2(l  ♦  +  —-M4  =  0  (a) 

n  n  GA  4  n  GA  n 


(4.8) 


(b) 


(4.9) 
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In  the  case  of  small  loss  tangents  of  order  (2.6)  it  can  be  shown  by 


straightforward  calculations  that: 


CRT 
R2  V 

C  “  PA 


(4.10) 


(4.11) 


<  <  ^p/^r2kEj/G»)t™VCcRV/<VkrlEK)  tMSE 


2fi£  kr2p/Gj-[l+an2Y2(l+kE^/G^)] 


It  is  seen  that  (4.10)  is  the  frequency  equation  of  an  elastic 
Timoshenko  beam  in  terms  of  real  parts  of  complex  moduli.  Its  validity  is 
based  on  the  usual  additional  assumption  that  real  parts  of  complex  moduli 

D 

vary  sufficiently  slowly  with  frequency.  Once  0>n  has  been  computed  from 
(4.10),  can  be  computed  from  (4.11). 

D 

A  mostly  sufficient  accurate  approximation  for  is: 


^  cVtx  -  Kv(i  ♦ 

gr 


(4.12) 


For  slender  beams  and  low  modes  the  first  term  in  each  of  numerator 
and  denominator  of  (4.11)  is  insignificant  relative  to  the  others.  Thus: 
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(a  2r2kER/GR)tan6  r  +  (cR  a_4/flR  )tan6„ 
n _ a  a  b  n  n  c 

1  +  on2r2  (1  +  kER/GR) 


(4.13) 


Substitution  of  (4.9)  into  (4.7)  results  in: 


-fl*t  xflRt 
wn(x,t)  -  Ansin  e  n  e  n 


(4.14) 


where  the  attenuation  is  now: 


a 

'n  n 


(4.15) 


To  obtain  an  idea  of  the  relative  importance  of  damping  due  to  shear 
and  rotator/  inertia,  the  attenuations  (4.5)  and  (4.15)  have  been  compared 
for  the  following  case:  Beam  of  rectangular  section 


l  a  40.0" 


h  ■  2.0" 


Material:  Boron  fibers,  Epoxy  matrix 

v-  *  v«  .5 
r  m 

Ef«60xl0*psi»  EmR«.SxlO*psi,  GmR».185xl0*psi 


tpn6_m  *  tan$  ®  *  tan6  »  .05 
cum 


ER  «  30.25  x  105psi  G.R  «  .544  x  lO'ps* 


tan6„ 

m 

155 
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tanS-  «  tan6_ 

u  m 


p  ■  1.78  x  10_<*  £b(mass)/ins 

1  !  : 

It  has  been  assumed  for  simplicity  that  real  paits  of  complex  moduli  and 

i 

loss  tangents  are  frequency  independent. 

For  the  first  mode:  .  .  1 

u»R  *  1480  1/sec 

l  I  ' 

}  Bending  only 

i 

ru  *  .308  1/sec  i  1 

fi>iR  -  1380  1/sec  i  i 

}  Timoshenko  beam 

» 

fii  *  4.44  1/soc  •  ' 

i  i 

It  is  seen  that  shear  and  rotatory  inertia  have  a  very  small  effect  on  the 

frequency  but  increase  the  attenuation  by  a  factor  of  14.4.  This  example 

1 

shows  that  for  damping  of  viscoelastic1  fiber  reinforced  beams  shear  and  rota¬ 
tory  inertia  are  of  major  importance. 

•  ;  .  I  ; 

i  v 

4.2  Forced  torsional  vibrations  of 

'  i 

laminated  cylinder. 

i 

A  thin  walled  cylinder  which  is  laminated  through  its  thickness  isj built 
in  at  one  edge  and  is  subjected  to  a  sinusoidal  forcing  torque  at  its 

other  edge.  Each  lamina  is  uniaxial ly  reinforced  and  has  the  same  material 

1  * 

properties  with  respect  to  its  material  axes.  The  laminate  is  symmetric- 

I 

with  following  lamination  scheme;  ,  . 

reinforcement  in  generator  direction  (axia?)  -  volume  fraction  vQ 


's 

£ 


2 

:? 


p 

% 

i 


H 

* 


3 


si 


reinforcement  in  +  6  direction  -  volume  fraction  v. 


reinforcement  in  -  8  direction  -  volume  fraction  v 


v8  *  v-6 


(4.16) 


For  the  purpose  of  analysis  of  torsional  vibrations,  the  only  effective 
laminate  property  needed  is  the  effective  complex  shear  modulus  Gu.  It 

i 

follows  from  (3.11)  that: 


i 1  -*  foV”  (+b\~ 

Gl2  *  Cj2 12  *  ‘,Ci2J2tQ/h  ♦  ^Ci212t+g/h  ♦ 


^^Gi212t_g/h 


(4.17) 


where  tQ  is  the  sum  of  the  thickness  of  the  0*0  laminae,  t+8  and  t_0  -  the 
sums  of  the  thicknesses  of  the  4-0  and  >0  laminae,  respectively. 


to/h  *  vq 


(4.18) 


t+0/h  *  l.Q/h  *  v0 


and  from  the  last  of  (3.16) 


(4.19) 


Introduction  of  (4.18-19)  into  (4.17)  yields 


G12  «  (0)ChhV0  ♦  2(9)Ch%v0 


(4.20) 


By  the  last  of  (3.16)  and  from  (2.4),  (4.20)  assumes  the  form: 

i 

Gi2  *  Gav0  ♦  [■|•(n+K-2X+GT)sin220+2GAcos^20]v0  21j 


!esa^saife»3..^£Sfa. 


cu 


It  has  been  mentioned  before  (Section  2)  that  the  imaginary  parts  of 
n ,1  and  k  can  be  neglected  with  respect  to  the  imaginary  parts  of  GT  and 
for  the  usual  fiber  reinforced  material.  Consequently,  the  separation  of 
(4.21)  into  real  and  imaginary  parts  assumes  the  form: 

Gi?  -  GjJ vq  ♦  [y(nR+kR-2£R+GR) sin220  ♦  2GRcos220]vQ 

(4.22) 

gTI  *  G^vo  ♦  (^G^sin220  +  2G*cos220)vQ  (b) 


In  the  usual  fiber  reinforced  materials  generally 


|(nR+kR-2AR+GR)  >  2GR 


I 

T 


<  2Gj 


It  follows  that  the  shear  loss  tangent  of  the  laminate 


Gu 

tand*  -  -r=r  (4.23) 

Gi? 

has  a  maximum  for  0®O  and  decreases  monotonically  to  a  minimum  for  0*45°. 

On  the  other  hand  shear  strength  is  smallest  for  0«O  and  increases  monotoni¬ 
cally  to  a  maximum  for  0®45#.  Therefore,  in  design  for  maximum  damping  it 
is  necessary  to  choose  the  smallest  angle  0  which  complies  with  allowable  shear 
stress . 

Let  the  forcing  torque  at  the  edge  xi«Jl  be  represented  as 

M  ®  M  sinwt 
o 


iik'Vr^'ielSSSSJ^^SSS^SJrSSSwSlJsJt&SS^^Ji*'* 


and  let  the  amplitude  of  angle  of  twist  <p  at  x»«i.  be  written  Amp  ($)• 
By  standard  theory  of  viscoelastic  vibrations  with  small  less  tangents 
[2,3; 


Amp  (<J>) 


cMq  ysin2(2a)+sinh2(2B) 

=  „*R  ,  cos  (2a)  +  cosh  (26] 
Jb|2w 


(4.24) 


where 

.2  _  G& 


0  *  o>lK6 
2c 

p  -  density 

C  -  Torsional  constant  of  section 
J  -  Polar  moment  of  inertia  of  section 
G*J  -  equ.  (4.22a) 

6*  -  equ.  (4.23) 

Numerical  analysis  has  been  carried  out  for  a  laminate  composed  of  boron/ 
epoxy  laminae  with  both  constituent  volume  fractions  equal  to  .5.  Laminae 


fractional  volumes  in  laminate  are 


v  «  .6 
o 


v±e  *  *2 


0  -  22.50* 


Analysis  has  been  performed  in  following  stages: 

(a)  Experimental  results  for  epoxy  matrix  complex  shear  modulus  and  loss 
tangent  as  a  function  of  frequency  have  been  described  by  an  empirical 
formula. 

(b)  This  formula  together  with  elastic  properties  of  fibers  have  been  used 
to  compute  effective  complex  moduli  of  the  uniaxially  reinforced  laminae, 
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as  a  function  of  frequency.  For  this  purpose  results  (2.7-9)  and  other 
formulae  given  in  [2,3]  have  been  used. 

(c)  With  the  aid  of  single  laminae  properties  the  real  and  imaginary  parts 
(4.22)  of  the  effective  complex  shear  modulus  Gf2  have  been  computed  as 
function  of  frequency.  It  should  be  noted  that  in  the  present  applica¬ 
tion  1  indicates  generator  direction  of  cylinder,  and  2  the  direction 
normal  to  generator  and  tangent  to  section  contour. 

*D  *  T 

(d)  The  results  for  Giz  and  G12  have  been  used  to  compute  (4,24)  as  a  func¬ 
tion  of  frequency  w. 

A  plot  of  such  results  is  shown  in  Fig.  3  for  a  cylinder  of  length  &=100  in. 
and  thin  walled  circular  section.  It  is  seen  that  the  first  resonance  peak  is 
very  significa.it  and  may  be  regarded  as  an  elastic  resonance.  However,  the 
damping  of  the  viscoelastic  matrix  becomes  more  effective  with  higher  order 
resonances,  the  fourth  one  being  considerably  reduced. 

5.  Conclusion 

It  has  been  shown  that  complex  moduli  of  uniaxially  fiber  reinforced 
materials  and  of  laminates  of  such  materials,  consisting  of  viscoelastic  matrix 
and  elastic  fibers  can  be  computed  in  straight  forward  fashion.  The  results 
can  be  used  for  analysis  of  structural  vibrations  on  the  basis  of  available 
theory. 

Two  structural  examples  have  been  given  to  assess  the  significance  of 
vibration  damping. 

Many  more  other  interesting  applications  can  be  analyzed  by  the  theory 
which  has  been  presented. 
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FIG.  1  UNIAXIALLY  FIBER 
REINFORCED  MATERIAL 


FIG.  2  LAMINATE 


